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This published work discusses problems arising from diffusion 
interaction between materials and their claddings and the film of con- 
densate formed on their surfaces. Since condensates practically always 
form on heated materials, the commencement of interaction is observed 
during the process of formation of the film. In this connection, the 
investigation of the characteristics of the diffusion interaction of the film 
with the materials should be divided into two stages: interaction during 

the process of the formation of the film and the interaction during the 
process of extended isothermal heating. Detailed results of two stages 
are presented in this work. 

A characteristic of the first stage is the simultaneousness of the 
course of the two processes: condensation and diffusional transfer of 

the cladding material into a growing layer. The continuity of the con- 
densate growth, as well as the law governing the growth of the layer 
itself, gives rise to the appearance of certain specific characteristics 
in the character of diffusion interaction (Section I). 

The interaction of the film with the material during isothermal 
heating may occur at the expense of a decrease of the content of one 
of the components of the cladding material and not as a result of the 
interchange reaction. It is borne out that geometric correlations have 
a substantial influence on the character of diffusion process in such a 
system. This is shown by the example of contact diffusion of monocar- 
bide of uranium with molybdenum and tungsten (Section II). 

Besides these general questions concerning characteristics of 
diffusion interaction of deposited materials, investigation is also made 
of the correlation of the diffusional mobility of atoms with thermodynamic 
characteristics of materials. 

In conclusion, a study is made of data concerning^ the diffusion of 
carbon in zirconium carbide (Section III). 

I. CHARACTERISTICS OF DIFFUSION PROCESSES IN CONDENSATES, 

FORMED ON A HOT CLADDING 

by. A. A. Ba bad-Zakhryapina and L. Gert 

1. Introduction 

To create a well bonded layer of condensate with a cladding 
surface, it is common practice to utilize heating of the cladding during 
the process of condensation. During the course of the whole process, 
there occurs an inevitable diffusion interaction between the cladding 
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material and the condensate material, a characteristic of which is a 
diffusional penetration of the cladding material into the layer of conden- 
sate. On the other hand, the penetration of the condensate material 
into the cladding also, to a certain degree, has to be determined by the 
conditions of condensation. 

The diffusion interaction between the cladding and the condensate 
may lead to a change in the conditions of condensation, i. e. , a change 
in the character of dependence between the thickness of the layer and 
the time it takes to form, as well as have a substantial influence on the 
phase formation in the condensate layer. Regardless of the fact whether 
these phenomena are desirable or undesirable in the process of formation 
of the condensate layer, they have to be taken into account and a method, 
based on prior knowledge, has to be developed by which the phenomena 
can be evaluated quantitatively. 

2. Unilateral Interaction of the Cladding Material with the Condensate 
Layer During the Process of the Layer Formation 

For the solution of this problem, it is necessary to know the 
growth law of the condensate layer. Subsequently, two laws will be 
examined: linear, i. e. , the growth of the layer with a constant speed 

and parabolic - with a variable velocity. 

The simplest case, which is encountered in the investigation of the 
diffusion interaction of the cladding with the condensate, is the unilateral 
diffusion in the condensate layer. This case may be formulated in the 
following manner: on a cladding of material A at a temperature T, sub- 

stance B is placed according to the law x = f(t) (t-time); it is necessary 
to determine the concentration of substance A (i. e. , c on the surface 
of the layer of condensate. 

Here it will be assumed that the coefficient of diffusion A in B 
(D^g) is considerably greater than the coefficient of diffusion of B in A 

(D^^). For simplicity and ease of notation, let us consider that 

D — O and D = D. 

BA AB 

Boundary conditions for this problem have the form (Figure 1) on 
the interface cladding condensate. 

d£ /9c \ 

C A dT = - D (*r) c ■ C(X ’ l) x = 0 = Co 

A 
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on the surface of the layer 


c(x, t) r = c A , where ft = f(t). 
x = IL A 

From the general solution of the Fick equation for a unilateral source 1 
and above boundary conditions fora case when ft = Zb\^Dt, we obtain: 

C A = 1 + F(b) (1) 

where F(b) = N/n-b exp b 2 >£(b), >3/(b) - function of Gauss errors. 


C 



Figure 1. Diagram of the Distribution 

of Concentration in a Layer of Condensate 
in a Case of Unilateral Single Phase 
Diffusion 


Under linear law of growth of layer, when ft - kt 
c 0 

C A = 1 + F(k, t) ^ 

where 

F ( k. ti = k^r « P -^*( 73 g r ). 

The expression, determining c under linear law of layer growth, 

A 

is an approximation and may be utilized for calculations at low values of 
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k and t. An accurate solution of this problem may be obtained by the 
solution of the equation of material balance on the boundary of cladding 
condensate, where 


- D 


9c 

9x 


x=i 


= c A k, 


(3) 


bearing in mind that dd /dt = k. 

Let us perform a substitution of variable k = d . Such a substitution 
is possible provided that 

9c _ dc A 

at £ - dt 


We will prove that the latter equivalence is valid only in the case if 
d = kt. Differentiating Equation (3) according to d and comparing the 
obtained result with the Fick equation for point d, we obtain 


d 

C A dl 



As c a ^ ® according tothe problem, then (d/d£) (d£ /dt) = 0 which gives 
d = kt. 


Bearing in mind that c (d , t) 
fer to variable n k M , we obtain 


c A and - D (9c / 9^) = c A k, after trans- 

A A A 


c(x, t) = Co exp 


( 4 ) 


(4) 


Let us prove that expression (4) satisfies Fick’s equation and boundary 
conditions. For this purpose, it is expedient to go to another system of 
coordinates X = d- x (Figure 2). In the new system of coordinates, 
expression (4) will be written as 


c(X, 


t) 


= co exp 


/kX - k 2 t 
\ D 



(5) 


4 




Figure 2. The Connection Between Systems 
of Coordinates Showing the Dependence 
of the Concentration of Cladding Material 
in the Condensate Layer on the Growth Rate 
Under the Linear Law for the Growth of the 
Condensate 

By direct substitution, it is easy to be convinced that expression (5) 
satisfies both the Fick equation and boundary conditions. 

From Equation (4), it follows that the distribution of the concentra- 
tion in the layer of the condensate, growing according to linear law, 
is pseudo- stationary. The form of the concentration curve in this case 
is determined only by the relationship k/D. 

With the increase of the growth rate, the concentration of layers at 
any one given point, and consequently on the surface, decreases. The 
same conclusion is arrived at by analysis of approximated cor relation(2) . 

In any case of parabolic law of growth, c remains a constant 

independent of time while its value is determined only by parameter 
"b". The above obtained correlations (3) and (4) maybe utilized for 
determination of coefficients of diffusion. The most suitable for this 
purpose is correlation (4). Conducting a linear interpolation of the 
concentration curve, which is permissible for lower values of kx/D, 
we obtain c(x, t) = c 0 D/(D + kx). In this relationship, magnitudes of n k M 
and "x" are easily determined experimentally. Thus, for computation 
of D, it is only necessary to know the concentration at any one given 
point along the thickness of the layer. In the case where, as a result 
of diffusion, solid solutions are formed and if the concentration depen- 
dence of the lattice parameter is known, then the concentration on the 
surface of the layer can be determined by X-rays. 
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Now let us examine a case where, as a result of diffusion, there 
is a possibility of formation of two phases in the growing layer: one 

with a large content of substance A (phase I) and another with a lesser 

content (phase II). The boundary conditions of such a problem appear 
as follows (Figure 3): at H = 0 c = Co, at ft ~ £1, there is a constant range 

of the concentration ci, 2 - c 2>1 at 1 = Jl 2 c = c . 



Figure 3. Diagram of the Distribution in the 
Condensate Layer in the Case of Two Phase 
Unilateral Diffusion 

This problem may be solved in the same manner as the preceding 
one but with one difference, that the solution is actually conducted 
separately for ranges 0 ^ x < SL\ and H\ < x Hz (£i and Hz are areas of 
existence of phases I and II) and is realizable as one single solution by 
means of the material balance condition on the boundary of the separation 
of phases. It is assumed that the movement of the boundary between 
phases or growth of phase I occurs according to the parabolic law 
xi - 2bi\/Di t. For the parabolic law of the layer growth (x 2 = 2b 2 \ZD 2 t) , 
there may be obtained an equation connecting the value of the concen- 
tration with parameters bi and b 2 : 

c 2, l - z _ c 2, l e F(bib z ) _ 1 

Co-ci >2 Co - Cl, 2 F(biCp) [ F (b 1 b 2 j + 1 ] F(bi) 

This correlation permits the evaluation of the character of the 
dependence between bi and b 2 . In the context of the problem, b 2 may 
vary within the limits bi < b 2 < oo. At b 2 = biF(bib 2 ) = 0 and 
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— C — = F(bi), i. e. , on the growing surface a constant concentration 

ci, 2 

c 2 i is maintained. 


If b 2 — oo, then F(bib 2 ) — oo and c 2> x - Ci, 2 


Co - c 2 i 


11 F ( b ) F(biCp) • 

this is equivalent to diffusion into semi-infinite space 1 . It is easy to 
demonstrate that bi exists only in a certain definite range of values and 
decreases with the increase of b 2 , reaching into limits of values, corre- 
sponding to a case of jet diffusion into semi- infinite space. 


Under the linear law of growth of the layer in an analogous manner 
(such as that for a case under parabolic law), an approximated correla- 
tion may be obtained: 

c 2, i ~ c XjZ _ c 2> i . F(ktbj) _ l 

Co- 01,2 c 0 - c i , 2 F(bxCp) [F(ktbi) + 1 ] F(bi) 


where 


F(ktbl) = ( kt - b '\[r^y*p(T \Jh ) * f tVdt)' 


to - time, 


upon conclusion of which phase II was formed in the layer. This corre- 
lation, just as for the case of one phase diffusion, may be utilized with 
small k and t. Its analysis shows that with the diminution of k and t, 
there occurs an increase of bj. The application of a quasi- stationary 
method of computation for the linear growth law permits computation 
of the initial single phase diffusion. The solution in differential form 
for this case has the appearance as follows: 



This system may be solved by any one of the well known methods of 
approximated solutions. As a result, it will be found that S.\ = Cp(t) 
and c A = f(t). 
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3. Interaction of the Condensate with the Cladding During the Process 
of the Layer Formation 


As in the preceding paragraph, an assumption will be made that 
in the process of the layer formation of condensate, there occurs a 
unidirectional diffusion of the condensate material into the cladding. 


Let us examine a case of parabolic law growth of the condensate 
and single phase diffusion (Figure 4). 

At xi = 2bi \[D \ t and x 2 - 2bWD 2 t, we obtain 


'B 


-ci 


ci 



NTiT bj exp b? [*(bi) + tt(b 2 )] 
*/tT b 2 exp b I ^(bi) + ^(b 2 )| 


where c 1 is the concentration of B in the condensate material of com- 
position A B , 
m n 

c is the concentration of B on the surface of the condensate, 

B 

Ci is the lowest concentration of B in the forming phase, 

Co is the concentration of B on the separation boundary of the 
cladding and the layer of the condensate, 
x 2 is the coordinate of the surface of the growing layer of the 
condensate, 

xi is the coordinate of the concentration stagger c* -0 of the 
forming phase. 


Estimation of the dependence c and bi on b 2 shows that an increase 

B 

of b 2 leads to a decrease of bi and an increase in c . At b 2 — ► °o , the 

B 

resulting expressions transform themselves into expressions describing 
diffusion from semi- infinite space 1 . 

However, if the velocity of the input of the substance into the cladding 
is so insignificantly small that it is immediately absorbed, then the prob- 
lem becomes considerably more complex. A variant of such a problem, 
namely, absorption of preliminarily accumulated layer, is investigated 
by Malkovich 2 . 
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Figure 4. Diagram of the Distribution of Con- 
centration During Unilateral Single Phase 
Diffusion of the Condensate Material into the 
Cladding 


4. Characteristics of Phase Formation in Condensate Layers Formed 
on a Hot Substrate 


It has been shown above that under parabolic law of growth of 
the condensate, the concentration of the cladding material on the surface 
of the layer remains constant; however, in the case of linear law, it 
varies with the passage of time in the direction of diminutions. In 
Figure 5, both laws of growth are shown in general form. At point A, 
the tangential line to the parabola is parallel to the straight line, 
corresponding to the linear law growth of the layer. On sector OA, the 
rate of the layer growth is greater than in the case under linear law, 
while on sector AB it is smaller. 


Let us juxtapose the characteristics of the layer growth of condensate 
c^ formed under parabolic and linear laws. With the thickness of the 

layer ft = OFgrown according to linear law, the concentration is greater 
on the surface. With the thickness OF < £<; Ocp, the surface concentra- 
tion will be smaller than c , i. e. , the rate of the layer formation on 

A 

this sector is greater than it is necessary for the support of concentra- 
tion c . on the surface. The ultimate thickness (limit) at which there 
A 


shall occur a drop in concentration from Co to c .is H - OCD. In a general 
case this will occur at OF < ft < Ocp during time t^< t< t^. 


It has been noted above that the highest rate of phase growth in the 
condensate is assured with the formation of the layer according to the 
parabolic law. However, in the case of formation of the phase during 
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diffusion into semi-infinite space, the phase growth occurs also according 
to parabolic law but with a lesser speed. 



Figure 5. Diagram Illustrating Correlations 
Between Rates of Growth Under Parabolic 
and Linear Law at Equal Time 


In Figure 6, two parabolas are shown cor responding to the growth 
of the phase in layer (1), formed according to parabolic law, and during 
diffusion into semi- inf inite space (2). 

With the growth of the phase in the layer formed according to linear 
law, it is possible to select such a rate of growth that in time t, the thick- 
ness of the phase in this layer will be lesser according to curve 2 in 
growth but greater according to curve 1, i. e. , St 2 < S. < £ 1# 

Thus, under the linear growth law of the layer, the picture of the 
distribution of phases is more extended than in the case of diffusion 
into semi- infinite space. This characteristic of phase formation under 
the linear growth law permits the growing of phases with a narrow area 
of homogeneity up to a certain noticeable thickness. 
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Figure 6. Diagram, Illustrating Correlations 
Between the Thickness of the Phases in the 
Case of the Growth of the Phase in the Layer, 

Formed According to Parabolic Law (Z), 

Upon Diffusion Into Semi- infinite Space (1) 
and in a Layer Formed According to Linear 
Law 

N67-34376 

II. DIFFUSION OF CARBON FROM CARBIDE OF URANIUM 
INTO MOLYBDENUM AND TUNGSTEN 

by V. S. Eremeyeiev, A. S. Panov, B. F. Ushakov, 
and Ye. V. Fiveiskii 


1. Introduction 

In connection with the search for new high temperature con- 
struction materials during the last two decades, many works have 
appeared on the investigation of the characteristics of carbides of 
transmutable metals 3 ’ 4 ’ 5 . Specifically, attention was focussed on the 
problems of diffusional mobility of carbon in carbides of transition 
metals 6 ’ 7 ’ 8,9 as well as the problems of contact interaction of carbides 
with metals 10 ’ 11 ’ 12 ’ 13 . 

In research on contact interaction of carbides with metals having a 
high melting point, two cases can arise. In the first case, an exchange 
reaction occurs of the type: 

Me^, C + Me^ = Me^ C + Me^ . (a) 
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In the second case, the cause of the formation of the carbide layer 
on the metal in contact is the migration of atoms of carbide into the 
metal. This migration is realized either as a consequence of surplus 
content of carbon in comparison with stoichiometry or because of the 
fact that carbide passes to the lower limit of the homogeneity area. 

With this, the following reaction occurs: 

ZMe T C + XMe TT = ZMe T C + XMe TT C , (6) 

I y + x II I y II z 

where C is the lower limit of the area of homogeneity of carbide 

Me C 

I y + x 

In principle, the reaction of the first type may be accompanied by 
the reaction of the second type. However, cases may occur when, in 
spite of the fact that the reaction of the first type is thermodynamically 
impossible, the contact interaction does take place. Systems UC -Mo 
or UC - W may serve as an example of this. Probably only by means 
of the reaction of the second type is it possible to explain the contact 
interaction of stoichiometric monocarbide of uranium with molybdenum 
and tungsten 14 ' 15 . 

The relatively limited available data from literature on the subject 
of the kinetics of solid phase reactions in such systems are often contra- 
dictory. Therefore, the study of contact interaction in system Me^C - 

Men, where the reaction of the type Equation (6) takes place, presents 

an unusual interest. 

In the current work, equations corresponding to the reaction of the 
type Equation (6) are solved, and discussion is presented on the influence 
of the carbon content on the kinetics of the process. For the purpose of 
substantiating the computed values, a study is made of the contact inter- 
action of the carbide of uranium with molybdenum and tungsten at 1000° C. 

2. The Solution of the Diffusion Equation for a Case of Migration of Carbon 
From Carbide into Metal 

The solution of a system of equations, describing the kinetics 
of solid phase reactions, often becomes complicated. Specifically, when 
the limiting stages are diffusion processes and the kinetic factors on the 
boundaries of separation of phases may be ignored, the finding of solutions 
is considerably facilitated. 

Let us examine the contact interaction in system Mc^ C - Me^. Let 
us limit ourselves by a case, when the kinetics of the growth of the 
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forming layers is determined by the migration of carbon from the carbide 
into the metal, and furthermore the reaction takes place according to the 
type Equation ( 6 ). The proposed scheme is applicable for those systems 
of Me^ C - Me^ when Me^ C is thermodynamically more stable than 

M e n C, while the growth of the carbide layer on Me^ is realized at the 

expense of excessive carbon in the carbide of Me^ C in comparison with 

the concentration of carbon corresponding to the lower limit of the area 
of homogenity of Me^ C. 

Figure 7 schematically illustrates the distribution of the concentra- 
tion of carbon in the process of jet reactive diffusion of carbon from 
carbide into the metal at moment of time t. 

In this connection, the following designations are introduced. C12 
and C22 are the concentrations of carbon on the boundaries of the form- 
ing carbide phase, C22 is the concentration of carbon corresponding to 
the lower limit of the area of homogeneity of the forming carbide, and 
C12 may be found after the solution. C33 corresponds to the maximum 
solubility of carbon in the metal and e is the thickness of the forming 
phase. 



Figure 7 . The Distribution of the Concen- 
tration of Carbon in the Process of Its 
Jet Diffusion from Carbide into the Con- 
tacted Metal at Moment of Time t 


Let us examine two cases; 

1 ) Diffusion of carbon from a semi- infinite body 

2 ) Diffusion from a finite body. 
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The first case belongs to the well known class of Stefan problems 
with parabolic law of growth 9 . 


a. The Case of the Semi- infinite Body 

For the three diverse phases Me^ C (phase 1 ), of the form- 
ing carbides Me^ C (phase Z) and the contacted metal (phase 3 ), let us 
put down the diffusion equation. 


ac. a 2 c . 

at ^i a x 2 
i = 1, 2, 3, 


(1) 


where index 1 corresponds to phase Me^ C, 

index Z is the forming carbide Me^ C, 
index 3 is the contacted metal, 

Di is the coefficient of diffusion of carbon in Me^ C, 

Dz is the coefficient of diffusion of carbon in Me^. C, 

D3 is the coefficient of diffusion of carbon in metal. 

Boundary and initial conditions have the form: 


Ci (x, 0) = C01 , x < 0 

Ci ( 0 , t) = Cn , 

C2 (e , t) = C22 > 

C 3 (€ , t) = C33 , 

c 3 (°°, t) = 0. 


( 2 ) 


The conditions of the mass balance on the separation boundaries 
at x = 0 and x = e are written down as follows: 



The separation boundary of phases at x = € is transposed according 
to the law 
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€ = 2a\/D 2 t . 


(4) 


Equation (4) determines the position of the phase boundary between 
the forming carbide and the metal in contact. The constant "a" is the 
means for the solution of Equation (1) with fulfillment of conditions (2) 
and (3). 


Omitting intervening insertions, we will set down the equation from 
which it is possible to determine "a": 



If D 3 < D 2 and C33 « C 22 , then the expression (5) may be simplified: 

V^F = * “p * 2 • (6 > 

Equalities (4) and ( 6 ) determine the speed of the formation of the 
carbide phase if the dimensions of the contacted carbide satisfy the 
conditions of the semi- infinite body. 

Let us note, however, that in case of fulfillment of equality ( 6 ), the 
concentration of C 12 in the forming carbide on the boundary of x = 0 is 
determined by equation; 


C12 = C22 ~ (C01 - erf (a). ( 7 ) 

a. 

If it turns out that C 12 is higher than the concentration of carbon C 12 , 
corresponding to the upper limit of the area of homogeneity in the form- 
ing carbide, then in order to determine the parameter of M a n instead 
of Equation ( 6 ), the following expression has to be used: 

* 

Cl2 ” C22 J r 2 /Q \ 

— = a v it erf a exp a . ( 8 ) 

C22 

Equality ( 8 ) follows from the assumption that on the boundary 
Me^ C - Me^ C, the maximum possible concentration of C 12 is maintained 
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while the kinetics of the growth of the layer are subject to parabolic 
law. 


The issue of formula (8) is given in work by Gertsriken and Dexter 16 , 
b. The Case of the Finite Body 


Let the diffusion of carbon occur from a plate of £ thickness. 
The initial concentration of carbon in the plate is equal to Coi. In the 
process of diffusion on the boundary of the plate at x = 0 and x = £, a 
concentration of carbon equal to Cn is maintained. Then utilizing the 
known solution for this case 17 , it is possible to demonstrate that the flow 
of the diffusion component on the boundary x = 0 and x - H, is equal to: 


P = 


4 Di 

i 


(Coi 


cii) E 

m = 0 


exp - 


Di (2m + l) 2 tt Z t 

r 


For simplification of the problem, let us examine a case when the 
flow of carbon, diffusing from the forming carbide into the metal, is 
insignificant and has practically no influence on the speed of the forma- 
tion of the carbide phase of the metal. 


We will find the speed of the growth of the carbide phase on the 
metal if on the boundary of the example at x = 0, a flow is unleased, 
determined by Equation (7) while at x = e the concentration is main- 
tained constant and is equal to Czz. In a general form, this problem 
is solved in other work . Utilizing this solution as applicable to our 
case, we may set down: 


y> 1 

k ^ 0 D ?(2k + ])| 


d 

d t< k + D 


2k + 1 



4 Di (Coi - 
a Czz 



exp 


Di (2m + l) 2 tt 2 t 

V 


Here the designations are the same, as those in point (1). 

First let us find the solution (8) in the assumption that in the right 
side we may be limited to one first member, corresponding to m = 0. 
This is assumed valid at 

Di tt 2 t _ , 
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The solution we seek is in the form of: 


e - Aq - Bq 


exp - 


Dx tt 2 t 

~T Z 


( 9 ) 


After transformation, the left side in Equation (8), taking into con- 
sideration Equation (9), is written as follows: 


s -2nat Y'' (2k + I). ■ ■ (2k - 2n + 3)Ap k ~ 2n + 2 Bo 2n * 1 [ (2n . k + 1 


n=: 1,2,... k = n, n+ 1 , . . 

00 00 

- 2 ( 2 n - l)a t V 


k + 1 


n = 1,2,... 


D 2 (2k + 1)| (2n)| (-1) 


(2k + i)... (2k - 2n + 3)Ao k ' 2n + 2 B 2 n ’ 1 


(2n - 1) 


al k + 1 


k = n - 1 , n, . . . 


D? (2k + 1)| (2n - 1) | (- 1 ) k 


U D 1 Tt* 

where a = z • 

a 2 

After substitution of sensible values £=0.4 cm, Dj ~ D 2 ~ 10" 10 
cm 2 /sec, Ao ~ Bo ~ 0.01 - 0.03, we obtain the information that the 
members of the row recede and disappear very quickly, and therefore 
it is possible with great accuracy to get by only with the first member 
(n = 1, k = 0), which is equivalent to: 

Bo Di TT 2 Di TT 2 t 

r exp 5 . 

£ 2 £ 2 

From this we find: 


Bo 


4 £ (C 01 - Cn) 

IT 2 C 22 


(ID 


For determination of 


= fJCoi " 
t = 00 2 C 22 

Finally, we set down 

, JMC 01 - C u ) 
2C 22 


Ao, we will make use of the obvious equality 

Cn) 


a 


TT 


exp - 


Dl TT 2 t 

£ 2 


( 12 ) 


(13) 
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The obtained solutions (4) and (13) complement one another. Solution 
(4) is valid at low temperatures and small "t M . Solution (13) is valid at 
high temperatures and high n t M . It is considerably more suitable to con- 
duct the investigation of the conditions of applicability of one form or 
another not in general, but with a specific system. 

Let us apply the above expressions for determination of layers of 
interaction in systems UC - Mo and UC - W. To bring the problem to 
the stage of designated layers, it is necessary to know the coefficients 
of diffusion of carbon in the carbides of uranium, molybdenum, and 
tungsten as well as in molybdenum and tungsten. 

3. Review of Data on Diffusion Constants in Molybdenum, Tungsten, 
and Carbides Mo 2 C, W 2 C, and UC 

Data on diffusion mobility of carbon in carbides of uranium, 
molybdenum, tungsten as well as in tungsten are collected in Table I. 

The coefficient of diffusion of carbon in the carbide of uranium at 
1000° C is obtained by means of extrapolation of data given by Samsonov, 
Strashinskaya, and Shiller 10 . 

The admissibility of such extrapolation is substantiated by experi- 
mental results 20 carried out in the temperature range interval 800° - 
1100°C. 

Literature contains no reliable data on the study of diffusional 
mobility of carbon in Mo 2 C and Mo at 1000° - 1400°. 

The results of work by Samsonov and Epik 7 concern much higher 
temperatures. 

Extrapolation of these data for 1000° C may lead to great error. 

Work by Takenchi, Homma, and Sato 20 describes the measurement 
of diffusional mobility of carbon in molybdenum with application of 
isotope C 14 . The authors of this work noted as well the fact of the 
formation of carbide of molybdenum during diffusional holding time, 
which lessens significantly the value of their results. The study of the 
diffusional mobility of carbon in molybdenum presents great difficulties 
in view of its low solubility in molybdenum. 

The experiments conducted by us according to the methods described 
by Gel'D and Lubimov 6 have shown that at 1000° C the coefficient of dif- 
fusion of carbon in molybdenum is no greater than 10 12 cm 2 /sec. Also, 
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taking into consideration the insignificant solubility of carbon in molyb- 
denum (about 0.02 weight percent C) zo , expression (5) for computation 
of "a" may be replaced by Equation (6). The same applies to tungsten. 

For determination of coefficients of the diffusion of carbon in M02C, 
the following work has been performed. Samples were prepared from 
melted down molybdenum (99. 9 percent pure) and placed in contact at 
1000° C with high purity graphite. 

The thickness of the M02C layer formed on molybdenum was deter- 
mined on a micrometic-durometer PMT- 3 as an average for 100 mea- 
surements. The results of the measurements are compiled in Table II. 

A constant concentration of C12 was maintained on the surface of 
molybdenum samples used in the conducted experiments. For this 
reason, Equations (4) and (8) may be used for computation of the 
coefficient of diffusion. The method employed for such computation is 
given by Gertsriken and Dexter 16 . The average value of the coefficient 
of diffusion of carbon in Mo 2 C after adjusting the experimental data 
(Table II) turned out to be equal to 2 • 10” 11 cm 2 /sec at 1000°. During 
computation 5 , it was accepted that C12 = 0* 536 g/cm 3 and C22 = 0.491 
g/cm 3 . 


4. Computation of the Thickness of Carbide Layers Formed on Molybdenum 
and Tungsten on Contact with Carbide of Uranium 


For computation of the thickness of forming carbide layers on 
molybdenum and tungsten after contact with carbide of uranium, for- 
mula (4) is utilized. Coefficient n a M is found from Equation (6). Equa- 
tion (6), as has been shown above, is valid, provided the flow of carbon 


under conditions (3.2) D3 



can be ignored. 


In our case, due to the low solubility of carbon in molybdenum and 
tungsten, this condition is fulfilled of which one may find substantiation 
using the following initial data: 


Cn = 0.6188 g/cm 3 at 1000°C 12 
C 2 2 = 0.4914 g/cm 3 for Mo z C 5 
C 22 = 0. 4671 g/cm 3 for W 2 C 5 
C 23 = 0. 001 g/cm 3 for Mo 22 . 
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Table I. Diffusion Constants in the Case of Diffusion of Carbon into UC, 
Mo 2 C, W 2 C, and W 


Diffusion 

Medium 

Method 

Temperature 
Interval, °C 

D 0 cm 2 / sec 

kcal/ mol 

D at 

0 

0 

0 

0 

0 



Authors 

UC 

c 14 

1470 - 2270 

2 • 1(T 2 

50 

5. 92 

• 10' 11 

9 

M 02 C 

M 

1400 - 1800 

780 

83 

5. 42 

• 10- 12 

7 

W 2 C 

M 

1525 - 1850 

25,000 

112 

1. 97 

• io - 15 

8 

w 

c 14 

1230 - 2014 

1.6 • 10" 6 ' 

1 

50 

4. 74 

• IO ' 15 

19 

Mo 2 C 

M 

1000 - 1400 

1 


2. 0 • 

io - 11 

Published 
W orks 


Notes: 1) C 14 - Method with application of radioactive isotope C 14 . 

2) M - Metallographic determination of the coefficient of diffusion along the 
thickness of the forming layer. 


Table II. Thickness in Microns of the Mo 2 C Layer 

Forming on Molybdenum After Contact with Graphite 


Temperature 

(°C) 

Time 

(hr) 

100 

250 

500 

1000 

8 ±2 

19 ±2 

28 ±3 


To elucidate the influence of carbon content in the monocarbide of 
uranium on growth rate of forming layers, let us select three values of 
concentration of carbon Coi in monocarbide equal to 0.624 g/cm 3 , 0.650 
g/cm 3 , and 0. 676 g/cm 3 which corresponds to 4. 8, 5. 0, and 5. 2 of 
weight in percent. 

We assume the density of the carbide of uranium to be equal to 
13 g/cm 3 . The values of the coefficients of diffusion are cited in Table I. 

The computed values of unknown constant "a" by means of the solu- 
tion of the transcendental Equation (6) are cited in Table III. 

The estimate of the concentration of carbon in the carbide of molyb- 
denum a| x = 0 according to Equation (7) shows that Ci 2 < Cu at 1000°C 
where Ci 2 corresponds to the upper limit of the homogeneity range of 
Mo 2 C. Therefore, it is necessary to utilize expressions (4) and (6) in 
the computation of the forming layer during the initial period. In the 
case of tungsten, it turned out that Cj 2 is greater than C* 2 and equal to 


20 





0. 55 g/cm 3 according to Samsonov 5 . Therefore, for tungsten, constant 
"a" is found from Equation (8). It is uniform for all compositions of 
the carbide of uranium 0. 3. 


Table III. Computed Values of Constant "a" 




Carbon Content in 



Monocarbide of Uranium, 


Temperature 

(°C) 

Percent of Weight 

Diffusion Pairs 

5.2 

5. 0 

4. 8 

UC-Mo 

1000 

0. 112 

0. 0612 

0. 0101 

UC-W 

1000 

1. 452 

1. 277 

0. 683 


Let us cite the final form of dependence of thicknesses of formed 
layers on time: 

€ = 9. 1 • I0~ 4f >jt forUC(4. 8 percent C) - Mo 

€ = 5. 5 • lCT 3 N/t for UC(5. 0 percent C) - Mo 

€ = 1.0* for UC(5. 2 percent C) - Mo 

€ = 2.7* I0~ 4 sft for UC(4. 8 - 5. 2 percent C) - W 

t in sec, € in p,. 

The solution of the equation of diffusion for a body of finite dimen- 
sions (Equation (13)), as has been noted above, is valid for extended 
times and high temperatures as well as with small thicknesses of 
In the investigated case at ft - 0. 4 cm, the parabolic law of growth should 
be replaced by the law determined by formula (13) after 10 5 hours. In 
this manner at a temperature of 1000° C in the wide time interval, the 
finite dimensions of the body may be ignored and the parabolic law of 
growth utilized. But, however, with high temperatures and thicknesses 
of samples about 0. 2 - 0. 1 cm, it is necessary to use formula (13), 
i. e. , to take into consideration the dimensions of the body already after 
10 - 100 hours. 

5. Comparison of Computed and Experimentally Obtained Values of Thicknesses 
of Formed Carbide Layers on Molybdenum and Tungsten 

With the aim of checking out the above presented computations, 
a task was set afoot to study experimentally the contact interaction of 
tungsten and molybdenum with the monocarbide of uranium. Samples 
of molybdenum and tungsten 09 mm and a height of 2 mm were cut from 
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forged rods obtained by the method of powder metallurgy. The purity 
of the metals was 99. 9 percent. The carbide samples of the two sub- 
stances (Table IV) 09 mm and h 4 mm were obtained by the powdered 
metallurgy metal. The density of the carbide samples was 12. 6 - 
12. 7 g/cm 3 . 


Table IV. Chemical Composition of Uranium Carbide 
(Weight Percent) 


No. 

of 

Parties 

U 

C Total 

C Free 

N 

O 

1 

95 

4. 7 - 4. 9 

0. 1 

0. 005 

0. 15 

2 

94. 8 

5. 0 - 5. 2 

0. 1 

0. 005 

0. 15 


The carbide and metallic samples were polished apd placed sequen- 
tially into molybdenum containers. 

The grooved connection of the cover of the container with the con- 
tainer assured a reliable contact. The diffusion annealing was conducted 
in a vacuum stove with a molybdenum heater in a vacuum 4 • 10" 4 mm 
of mercury column at 1000°C up to 2500 hours. The accuracy of the 
temperature measurement was ±20°C. After annealing, the samples 
from tungsten and molybdenum were investigated rontgenographically, 
metallographically, and radiometrically by measurement of a- activity. 
Diffusional pairs were used for computation because they possessed 
satisfactory linkage. The experiments were duplicated. 

The penetration of uranium into the contacted samples was insigni- 
ficant for all investigated compositions. Along the layer, radiometric 
analysis did not uncover traces of uranium at a depth of 20 p after 
diffusional annealing during 2500 hours. Therefore, it may be considered 
that the processes of mass transfer in the first approximation are deter- 
mined by the diffusion of carbon. 

On the surface of metallic samples after diffusional annealing, the 
layers of the resulting interaction were determined and measured 
metallographically. 

The thickness of the layers is presented in Table V. 
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0 500 1000 1500 2000 2500* t(hour) 


Figure 8. The Computated and Experimentally Obtained 
Values of the Thickness of the Formed Carbide Layers 
on Molybdenum After Its Contact with Carbide of Uran- 
ium at 1000° C. Curve 1, Obtained by Computation, 
Corresponds 5. 2 Weight Percent C; Curve 2, 5. 0 Weight 
Percent C; and Curve 3, 4. 8 Weight Percent C in Carbide 
of Uranium. 


Author's results: 5. 0 -f 5. 2 Weight Percent C in UC 

4. 7 f 4. 9 Weight Percent C in UC 

Data from literature: x 5. 1 Weight Percent C in UC 9 



Figure 9. Computed and Experimentally Found Values 
of Thicknesses of Formed Carbide Layers on Tungsten 
After Its Contact with Carbide of Uranium at 1000° C. 
Curve 1 Corresponds 4. 8 - 5. 2 Weight Percent in UC. 

Author's results: 5. 0 - 5.2 Weight Percent C in UC. 
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Table V. Results of Interaction of Molybdenum and Tungsten 
with Carbide of Uranium at 1000° C 



Time 

Thickness of Interaction Microns 

No. of Parties 

(hr) 

Mo 

W 

1 

(4. 7 - 4. 9 wt %) 

1300 

2 - 3 

0. 5 

2 

(5.0 - 5. 2 wt %) 

250 

5 

0. 5 


WEM 

8 

0. 5 


mSm 

10 

0. 5 


EH 

+ 17 

1 4 2 


Rontgen analysis has shown that the formed layers correspond to 
M 02 C for molybdenum and WC, W 2 C for tungsten. In tungsten, in all 
cases X-ray analysis has shown the presence of phases WC and W 2 C. 

As long as the sensitivity of the metallurgical method is not greater 
than 0.5 -r 1 micron, (then) the layers of carbides of tungsten were not 
noticed after extended contact of 250, 500, and 1000 hours with UC 
(5. 0 - 5. 2 weight percent) and 1300 hours with UC(4. 7 - 4. 9 weight 
percent). 

Experimental data are presented in Table V together with computed 
values in Figures 8 and 9. In these very diagrams are also presented 
the results of work by Katz 11 . The computed and experimental values 
are in excellent agreement. Therefore, the proposed diffusional dia- 
gram of contact interaction is fulfilled, at least at 1000° C. 

6. Conclusions 

1) An equation for diffusion is solved which corresponds to 
the migration of carbon from carbide into (he metal in contact with sub- 
sequent formation of carbide layer. 

2) The interaction of carbide of uranium with molybdenum 
and tungsten at 1000° C up to 2500 hours has been experimentally 
investigated. 

3) The computed values of the magnitudes of the forming layers 
are in satisfactory agreement with those arrived at experimentally. 
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III. RESEARCH INTO THE DIFFUSION OF CARBON INTO ZIRCONIUM 
AND ZIRCONIUM CARBIDE 

by R. A. Andrievskii, V. N. Zagryazkin, 
and G. Ya. Meshcheryakov 

Information on the diffusional mobility of carbon in zirconium and 
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its carbide is very limited. Babikov and Gruzin , studying the electro- 
transfer of carbon into zirconium at 900°- 1200° C, have evaluated the 
numbers (numerical values) of the transfer. The determination of 
effective charge was not possible due to the absence of data on diffusion. 
In the work of Samsonov and his colleagues 7 , the diffusional saturation 
of zirconium by carbon was investigated, and certain energy parameters 
of interaction were evaluated. As far as it is known, these works force 
the only information on the diffusional mobility of carbon in zirconium 
and carbide of zirconium. 

The samples of iodide zirconium and carbide of zirconium used by 
us in our work had the dimensions 10-12 mm and h 8 - 10 mm. 

The latter sample was prepared by means of hot pressing at 2200°C 
under pressure of about 200 kg/cm 2 with a subsequent two hour anneal- 
ing under the same temperature for homonization. The porousness of 
the samples did not exceed 5 percent; the composition of bonded carbon 
was 11.2 percent, which corresponds to the formula ZrCo. 96! the quan- 
tity of free carbon was 0.65 percent, zirconium - 88.0; the lattice para- 
meter - 4. 697 A. The iodide zirconium contained v weight percent) 

0.01 C; 0. 035 N; 0.11 O; 0. 22 Fe; and 0.015 Hf. he quantity of the 
other admixtures was less than 0.005 weight perci it. 

The radioactive carbon C 14 was smeared in a thin layer on thoroughly 
polished bars of samples in the form of an emulsicn in glycerine 
(specific gravity 30 microcurie/g). The samples vere annealed by pairs 
in a graphite container which assured a tight contact of "active surfaces". 
The zirconium was annealed in a vacuum 10” - 10 5 Torr at 1100° - 

1600° C and the carbide of zirconium was annealed in refined helium at 
1600° - 2100° C. The temperature was measured by thermocouples 
Pt + 10 percent Rh/Pt and W + 5 percent Re/W + 20 percent Re with an 
accuracy ±20°C. After annealing, the samples were separated and a 
layer radiometric analysis was performed with construction of depen- 
dency lgj^= f(x 2 ), where ^ is the intensity, p is the numbering (score), 
and x is the distance. In all cases, these dependencies had a linear character 
which evidenced adherence to conditions of diffusion from a thin layer 
(Figure 10). However, for zirconium (the initial one), two points some- 
times do not fall on linear dependencies lg/-x 2 , but microscopic and 
microdurometric investigations did not show formation of a carbide 
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phase. The coefficients of diffusion D were computed according to a 
well known formula: 

^ 0.1086 

D = — 

r tan a 

where r is the annealing time and tan a is the tangent of the angle of 
incline of the straight (direct) line in coordinates lgY-x 2 . The method 
of the least squares was used for calculation of data. 

Tables VI and VII and Figure 11 show the values of the coefficients 
of diffusion of carbon in zirconium and the carbide of zirconium. All 
values are averages from 3-4 measurements. For the diffusion of 
carbon in the carbide of zirconium, another method was also used. It 
consisted of the application of the carbide of zirconium containing radio- 
active carbon- 14. The thickness of such a sample was 6-8 mm. 
Computation was conducted according to the well known formula for a 
case of diffusion from a constant source into infinite space. Both methods 
lead to very similar results. The difference in values of D fell well into 
the margin of error of the experiment. The maximal error in the deter- 
mination of the coefficient of diffusion and energy of activation in our 
case consisted of about 20 and 4 percent, respectively. 


a b 



Figure 10. Results of Layer Radiometric Analysis in Coordi- 
nates ig/ -x 2 : a - Zirconium; b - Carbide of Zirconium 
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The temperature dependencies of D, obtained with the help of the 
method of least squares, has the form 

D =3.6* 10' 2 • exp (-34,200) RT/cm 2 /sec (1) 

c * r 

D „ „ =3.3- 10" 2 • exp (-114,000) RT/cm 2 /sec. (2) 

,c— *-ZrCo.96 

Let us discuss the obtained dependencies (1) and (2). First of all 
it is expedient to evaluate the entropy of activation AS of diffusion of 
carbon in (3-Zr, utilizing the well known expression: 



where Do is the preexponential factor; p is the number of equivalential 
transpositions of interstitial atoms, equal to 4 for b. c. c, (body centered 
cubie); v is the frequency of fluctuation, equal according to evaluation 
(Q/2m X. 2 ) 1 / 2 ; Q is the energy of activation; m is the mass of interstitial 
atoms; \ is the distance between interstitial atoms; a = 1/24; a is the 
grid constant; and R is gas constant (3). For values Do = 3. 6 • 10“ 2 
cm 2 /sec, v = 1.35 * 10 13 sec -1 , and a = 3. 616 A, we obtain the entropy 
of activation, equal to 4. 95 cal/g-mol • °C. It was not possible to 
compare this value with the computed one because data on temperature 
variation of the Young Modulus for (3 of zirconium is not known. 


Table VI. Coefficients of Diffusion of Carbon 
in (3 of Zirconium 


Conditions of Experiment 


Temperature 

(°C) 

Duration 

(hr) 

D * 10 7 , cm 2 /sec 

1100 

2 

1.4 

1200 

5 

2. 2 

1300 

1.5 

7. 3 

1400 

4 

15 

1520 

2 

27 

1600 

0. 5 

33 
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Table VII. Coefficients of Diffusion of Carbon 
in Carbide of Zirconium 


Conditions of Experiment 


Temperature 

(°C) 

Duration 

(hr) 

D , c m 2 /sec 

1580 

7 

9. 5 • 10 -12 

1650 

3. 5 

rH 

r-H 

O 

00 

CO 

1780 

3 

4. 1 * 10- 10 

1900 

2 

2 • 10 -9 

2080 

i 

2 

i 

6 • 10- 9 



4 5 6 7 


I0 4 /T°K 

Figure 11. Temperature Dependence of Coefficients 
of Diffusion of Carbon in Zirconium (Equation (1)) 
and Carbide of Zirconium (Equation (2)) 


It is expedient to juxtapose the parameters of diffusion of carbon 
in (3 of zirconium obtained by us with analogous data for other elements 
of penetration ( inter stition) . Table VIII presents these data. 

Bearing in mind the correlation of Wert and Zener 24 , 


AS = - Q 

3T 


(4) 


28 





where \x and p 0 are modulus change at a given temperature and 0°K, it 
is evident that the relation of entropy of activation to energy of activation 
should be constant for the diffusion of various interstitial elements into 
the same metal. For convenience, this relation (AS/Q) is usually multi- 
plied by and a dimensionless parameter is obtained 25 , characteriz- 

ing the temperature variation of the shift module. 


AS 8(n/Ho) 

Q PL 3 ( T / T ) 

F .L 


(5) 


Expressions (4) and (5) assume the absence of relaxation which, of 
course, is only the first approximation. 

Table VIII presents values of AS and (3 computed by us according to 
other data 26 ’ 27 ’ 28 ’ 29 . As it is apparent from presented results, values (3 
for diffusion of hydrogen, nitrogen, oxygen, and carbon in zirconium 
have basically the same order. The concurrence of results for oxygen, 
nitrogen, and carbon are relatively good, but somewhat poorer for 
hydrogen. In absolute magnitude, the values evaluated by us had a 
meaningful character. Thus for other metals, the b. c.c. parameter 
p, determined from direct measurements of temperature dependence 
of the module change, comprises, for example, 0.43 for a iron and 
0. 40 for tantalum 30 . 


Analysis of the results of Table VIII thus shows that the parameters 
of diffusion of various interstitial elements into p of zirconium are 
juxtapositional. 

In the investigation of diffusion of interstitial elements, it is common 
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practice to use a model of "touching hard spheres” * ’ . In work by 

Schumann , attention was directed to the fact that the energy of the 
activation of diffusion of hydrogen, nitrogen, and carbon in a iron is 
dependent linearly on the radii of these two elements, and, furthermore, 
the "zero” energy of activation corresponded to the dimension of the 
octahedral pore. 

An analogous dependence was observed as well for diffusion of 
interstitial elements into p of zirconium (Figure 1Z). The atomic radii 
r^. = 0. 36, ro = 0. 6, r^ = 0. 70, and r = 0. 77 A are taken from other 

works 13, 14,15 and relate to covalent radii. 

It is characteristic that the value of the radius, corresponding to 
the "zero” value of activation energy, comprises about 0.25 A which is 
very similar with the radius of the octahedral pores in p of zirconium 
( r OKT = 24 A). 
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The first direct indications on the location of interstitial elements 
in the [3-zirconium lattice are absent. As it is known, for other b. c. c. 
metals, the body centered cubic based on the results of internal friction 
and electromicroscopic investigations, the interstitial atoms are located 
primarily in octahedral holes J . 



20 40 60 80 


o 

A 

Figure 12. Dependence of Activation Energy of Interstitial 
Elements in (3 of Zirconium on Magnitude of Atomic Radius 


For hydrogen atoms, however, a more predominant distribution in 
tetrahedral holes seems more natural. This is possibly explained by a 
higher solubility of hydrogen in [3 phases of transition metals of group 
IV in comparison to hexagonal a phases. The dimensions of tetrahedral 
holes in the b. c. c. is considerably greater than in hexagonal. The 
estimates made by Magnier and Accary 12 also add evidence to the fact 
of predominant distribution of hydrogen atoms in tetrahedral holes of 
the b. c. c. , of metals. 

The circumstance, that regardless of possible distribution of 
interstitial atoms in octahedral (O) or tetrahedral (T) holes of [3- 
zirconium lattice dependence Q = f(r) has a linear character, signifies 
that the diffusion is connected with transition through octahedral holes, 
i. e. , it can be represented by diagrams 0-0 or T-O-T. From this it 
follows that transitions 0-0, as has been noted by Magnier arid Accary 12 , 
may be realized according to diagram O-T-O which, evidently, is 
characteristic for diffusion of carbon, nitrogen, and oxygen. 
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Absence of data on the modulus change for monocrystal (3 of zirconium 
does not permit the qualitative evaluation of the possibility of filling the 
tetrahedral and octahedral holes by interstitial atoms. 

However, conclusions based on examination of diffusion of inter stitial 
elements within the framework of a model of ’’hard spheres” are envisaged 
as not being unquestionable. 

In such an approach, the electronic structure of solid interstitial 
solutions is not taken into consideration. The real value of atomic radii 
in solid solutions still remains unknown. Evaluations carried out accord- 
ing to formulas of Pauling and E. S. Sarkinson 14 * 15 indicate that the radius 
of a hydrogen atom may increase in solid solutions up to approximately 
50 percent. The increase of atomic radii of other interstitial elements 
is considerably smaller (up to approximately 10-20 percent). Further 
investigations must refine the real character of dependence Q= f(r) 
and 1 ’ life span” of interstitial atoms in tetrahedral and octahedral holes. 

The theory of the diffusion of metallic and metalloid atoms in inter- 
stitial phases are not developed yet. Therefore, discussion of dependence 
(2) may be rather qualitative. Thus, it is not possible to explain the 
rather high significance of the preexponential factor. 

In contrast to the earlier examined case of diffusion in p of zirconium, 
the diffusion of carbon in the carbide of zirconium is realized not accord- 
ing to the principle of interstitial but according to the principle of sub- 
stitution in the metalloid sublattice of this combination. In this connec- 
tion, taking into consideration as well the difference in the melting 
temperatures, the considerably lower diffusional mobility of carbon in 
the carbide of zirconium and considerably higher energy of activation, 
in comparison with diffusion of carbon in (3 of zirconium, seem to be 
absolutely natural. 

1 8 

Borisov, Lubov, and Temkin made a study of the evaporation of 
the carbide of zirconium and established that around composition 
ZrCi.oo the evaporation occurs congruently. The value of the partial 
heat of sublimation Q of carbon, according to the data of this work, 
comprised about 178 kcal/mol. 

For high melting transition metals, the relation Q/Q is usually 
0. 6 t 0.8 19 . In our case, this relation turned out to be equal to 0. 65. 

Bearing in mind this circumstance as well as the high melting 
temperature of ZrC, the obtained value of the activation energy for 
carbon diffusion in zirconium carbide appears to be quite reasonable. 
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